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It is challenging to implement genuine free running single photon detectors for the 
1550 nm wavelength range with simultaneously high detection efficiency (DE), low 
dark noise, and good time resolution. We report a novel read out system for the 
signals from a negative feedback avalanche diode (NFAD) 1 3 which allows useful op- 
eration of these devices at a temperature of 193 K and results in very low dark 
counts 100 CPS), good time jitter (<~ 30 ps), and good DE (^10 %). We char- 
acterized two NFADs with a time correlation method using photons generated from 
weak coherent pulses (WCP) and photon pairs produced by spontaneous parametric 
down conversion (SPDC). The inferred detector efficiencies for both types of pho- 
ton sources agree with each other. The best noise equivalent power of the device 
is estimated to be 8.1 x 10 -18 W • Hz -1 / 2 , more than 10 times better than typical 
InP/InGaAs SPADs show in free running mode. The afterpulsing probability was 
found to be less than 0.1% per ns at the optimized operating point. In addition, we 
studied the performance of an entanglement-based quantum key distribution (QKD) 
using these detectors and develop a model for the quantum bit error rate (QBER) 
that incorporates the afterpulsing coefficients. We verified experimentally that using 
these NFADs it is feasible to implement QKD over 400 km of telecom fibre. Our 
NFAD photon detector system is very simple, and is well suited for single-photon 
applications where ultra-low noise and free-running operation is required, and some 
afterpulsing can be tolerated. 
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I. INTRODUCTION 



Single photon detectors (SPD) are important devices in many research fields, such as 
optical quantum computing^, quantum cryptograph}^, spectroscopy, fluorescence lifetime 
measurements 7 , and light detection and ranging (LIDAR) 8 . The telecom wavelength range 
around 1500 nm is of particular interest for various reasons. For example, quantum key 
distribution (QKD) via optical fiber must operate at the lowest photon transmission loss at 
1550 nm in order to achieve the longest distance. 

SPDs can be operated in two distinct modes: gatecP or free running^^l. Free running 
detectors are very important in certain applications, particularly where photons arrive in a 
random fashion, rather than synchronized on a regular clock time basis. For example, this 
is the case for photon pairs produced by spontaneous parametric down conversion (SPDC) 
with a continuous wave pump. 

In the telecom wavelength range there exist several good candidates for free running SPDs, 
including photomultiplier tubes (PMT)P^, passively quenched InP/InGaAs SPADs 16 , shallow 
junction silicon SPADs using frequency upconversion 17 , transition edge sensors (TES) 18 , and 
superconducting nanowire single photon detectors (SNSPDs) 19 . 

The SNSPDs currently have the best overall performance, such as low darkcount rates 
(DCR < 100 Hz) and good timing jitter (< 65 ps). However, owing to their superconducting 
temperature requirements (under 4.2 K), cryogenic systems are required. Moreover, the de- 
tection efficiency (DE) of SNSPDs is dependent on the polarization of the input photond^SED 
On the other hand, negative feedback avalanche diodes (NFADs), a recently developed type 
of InP/InGaAs APD (Princeton Lightwave Inc), can be operated in free running mode with- 
out cryogenic temperatures. These devices possess an integrated passive quenching resistor 
that minimize the amount of avalanche charges produced by a photon detection (another 
approach to self-quenching in InP-based single photon APDs, which employs epitaxial bar- 
riers as a negative feedback mechanism ! 22 1 23 1 ) . Furthermore, these APD detectors have no 
DE dependency on the input photon polarization. 

The darkcount rate and timing jitter of NFAD detectors arise from some intrinsic mate- 
rial properties and the device structure. Their darkcount rates are induced by the size of 
their avalanche region, as well as the trap-assisted tunneling effectP (TAT). Timing jitter is 
generally dominated by the stochastic nature of the avalanche breakdown process. However, 
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the external operation setup can be optimized to achieve a high signal-to-noise ratio (SNR) 
and good time resolution. We developed a novel readout system for the NFADs, which 
features (a) very low operating temperature; and (b) minimal avalanche charges required to 
trigger the readout circuit. This allows our system to achieve ultra low darkcount rates. 

The significant features of our readout method are: first, the NFAD operates in a genuine 
free running mode; there is no active quenching circuitry involved to impose additional 
deadtime (which is usually needed to suppress the afterpulsing effect and to lower darkcount 
rates); Second, a pulse transformer is employed to isolate the NFAD bias circuitry from the 
high speed RF amplifiers to read out the photon electric signals. Instead of conventional 
capacitive coupling, we use inductive coupling, so that only the variation of the avalanche 
current is coupled to the external amplifier. This dramatically improves the sensitivity of 
the readout. 

In this work two different NFAD devices are studiecP. Detector I/NFAD 1 (model no. 
E2G6) has a feedback resistor of 1.1 and a diameter of 22 fim. Detector II/NFAD 2 
(model no. E3G3) has the same nominal feedback resistor, but a larger device diameter of 
32 /im. 

The paper is organized as follows: we start with the experimental setup and device details. 
Then the "Simulation Program with Integrated Circuit Emphasis" (SPICE) model of our 
readout circuit including the NFAD equivalent circuit is simulated and compared with the 
experimental results. This is followed by a measurements of the DE, which is later verified 
using two NFADs to detect photon pairs generated by SPDC. Next the results from the 
afterpulsing, deadtime, and timing jitter measurements are presented. Lastly, we study the 
feasibility of using our NFADs to perform entangled photon QKD over distances of more 
than 400 km. 

II. THE DETECTOR SYSTEM 
A. Setup 

Fig. [T] depicts the detailed design of our in-house developed readout circuit. The entire 
system is comprised of two parts. The NFAD and its bias circuitry is inside a deep freezer 
(Global Cooling, ULT-25) with a lowest base temperature of 187 K. At this base temperature, 
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FIG. 1. The system design diagram of the entire NFAD readout circuits and connections. All the 
components in the right panel are inside a deep freezer. The ones in the left panel are at room 
temperature. SMA is the high frequency connector; FC/PC is the single mode fiber connector; 
Hirose is the NFAD high voltage cable connector. 

the detector operating temperature is 193 K. The primary coil of the pulse transformer (Pulse 
Electronics) is connected between the system ground and the cathode of the NFAD. The 
secondary coil of the transformer couples to a 50 Q transmission line (TL). A bypass resistor 
of 10 MQ and a Schottky diode serve to protect the NFAD from unexpected voltage surges 
or spikes. 

The rest of the readout electronics, including high voltage modules and the wide band 
RF amplifier, are in a nuclear instrumentation module (NIM) at room temperature. The 
high voltage control module (Matsusada) supplies variable bias voltage to the NFAD, and 
its output value is controlled by a multi-turn linear potentiometer. The bias voltage spans 
from to 80 V. Two wide band RF amplifiers (Mini-Circuits) have 40 dB gain in total to 
amplify the signals from the pulse transformer via the TL. 

The experimental setup to characterize the NFADs is based on time-correlated single 
photon counting (TCSPC), see Fig. [2j Optical pulses are generated with an id-301 1550 nm 
distributed feedback (DFB) laser, attenuated by an EXFO programable attenuator (FVA- 
3100). The laser is triggered by a two-channel function generator (Tektronix AFG3502). The 
attenuated weak coherent pulses pass through the attenuator and are sent to the NFAD. 
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FIG. 2. Testing setup for characterization of the NFAD using the time-correlated single photon 
counting method. Details in the text. 



The time tag unit records detections and the synchronization signals, and establishes their 
coincidences. 



B. The photon detection pulse responses 

The photon detection response is measured by a fast digital oscilloscope. Fig. [3] (a) shows 
a typical waveform from an experimental measurement of the single photon response. We 
also simulated the response using a SPICE model for the equivalent circuit depicted in Fig. 
[3] (b) and the NFAD model based on the results of Hayat et al. 24 . The simulation allows us 
to estimate the heating induced by each photon detection event, as well as the amount of 
avalanche charges coupled to the external readout circuit. Both estimations can be achieved 
by simulating transient response simulation of a single avalanche event. 

The simulation results for the signal response of two readout methods are illustrated in 
Fig. [3j (c) shows the transient response based on our readout circuit design (see Fig. [3] 
(b)); graph (d) is the result for the conventional readout method shown in Ref. 1. The 
conventional circuits transient response has a persistent current plateau ranging from ^10- 
100 ns and is proportional to the overbias voltage 25 . In comparison, our main response 
voltage peak is less than 5 ns, and is independent of overbias voltage. The slow variations of 
the persistent current induce almost no voltage via the transformer to the amplifier input. 
The simulation and experiment results agree well for the leading period of about 5 ns long. 
The experimental data exhibit a more complicated response, because the current SPICE 
model only accounts for the main physical process. 

In order to study the power consumption arising from the Joule heating effect, we per- 
formed numerical integration over the product of total power supply current and volt- 
age. Two sources of heat generation are identified: a constant background heat power 
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FIG. 3. (a) The experimental single photon response recorded with a fast digital oscilloscope 
after the 40 dB gain amplification, (b) The SPICE model for the NFAD read out circuit. The 
dash line square highlights the NFAD equivalent circuit components; The testing points TPvl and 
TPdvl represent photon detection response output and intrinsic SPAD voltage drop, respectively; 
the voltage source VBR represents the breakdown voltage of the NFAD. (c)-(d) are results of the 
transient simulation of a single photon responses, before the amplifier, (c) is the simulated response 
for the readout circuits used in this work; (d) is the simulated response with the conventional 
readout methods using a bias-tee in Ref. 1, The model parameters for the bias tee in the the 
conventional method are simplified from the ones provided by its manufacturer. 

of r>j 5 x 10 -4 W induced by bias and protection circuitry (e.g. 10 resistor Ri), and an 
additional 4 x 10 -12 J produced by every detection. Therefore, even at very high detection 
rates, the heating caused by detections is negligible, compared to other sources of heat, such 
as the background heat. 

Moreover, it is of interest to analyze the net avalanche charge, Q av , coupled to the exter- 
nal amplifier. We performed the numerical integration over the simulated voltage across a 
50 ft load via Qqm — Jo Vtpvi • dt. The computed result showed the avalanche charge with 
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inductive coupling is about 5 x 10 -16 C (corresponding to 3 x 10 3 electron charge e~). Com- 
paring this to the net avalanche charge for a conventional readout circuitP, ours consumes at 
least two orders of magnitude less avalanche charges to produce a legitimate photon response 
pulse. This low coupling charge leads to the improved sensitivity of our readout circuit, and 
contributes to the low afterpulsing dark noise performance of the detector system. 



III. DETECTION EFFICIENCY CHARACTERIZATION 

A. Repetitive weak coherent states method 

We first use attenuated coherent states, prepared at regular repetition rate, to characterize 
the detection efficiency of the NFAD. We choose the period of the repetitive pulse (T = 
100/xs) to be much longer than the detector recovery time (r = 50 ns). The coherent state 

| a |2 OO ^ 

of each pulse is represented in the Fock basis by \a) = e~^~ ^j== where the average 

71=0 

photon number contained in each pulse is \i = \a\ 2 . Since detection events are identified 
by comparing the voltage level of an photon response pulse to a pre-setting threshold, we 
model this type of detector as a so-called bucket detector 26 . The probability for click and no- 
click of the bucket detector are described by / C H C k = 1 ~~ (1 ~~ Vd)™ \n) (n\ and / no click = 

77=0 

OO 

^ (1 — r]E>) n \n) (n\ respectively. 

77=0 

In the experiment, we measure the ratio of detected photon click events, i?det, within 
the 1 ns correlated time window, r , versus the photon trigger repetition rate, / repj under a 
certain average photon number per pulse, \i. We combine the coherent quantum states to 
derive the detection efficiency, rju with the TCSPC method as: 

where the total non-photon detection events, g?b, becomes a summation of pure dark counts, 
L>, and afterpulsing events induced by photon detections, rf a f p , i.e. = D + rf a f p . r is 
the smallest time window in the TCSPC counting system. Here we always use a 1 ns time 
window for all measurements with the time tag unit. 

Because the darkcount and repetition rates are both below 10 kHz, and r = 10 -9 s, we 
are able to use 1 — ^b^o — 1 as a very good approximation in most of our DE measurements. 
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Thus, we have a practical relation to compute the DE: 

tto*~ln(l-^Y (2) 

V \ /rep / 

Here we assume the detector operator is known to behave as a bucket detector. If one wants 
to characterize a completely unknown detector, then a more general approach would be to 
perform the detector tomography^ZEEl 

Prior to the detection efficiency (DE) measurements, we first measure the darkcount rate 
for various bias voltages. Both detectors exhibit a similar behavior of voltage to darkcount 
rate. Fig. [| shows the trend of darkcount rate as a function of input voltage. Over the 
range of overbias voltages tested, the darkcount rate varies from 0.1 to 10,000 CPS. 

The DE testing conditions are as follows: the laser source, id-301, is triggered at 1 kHz 
with its output attenuated such that the photon flux is determined to be 1 photon per 
pulse by a calibrated power meter 29 . We then verify this attenuation using a commercial 
detector^, id-201, set to 10 % detection probability with a 20 ns gate and a 40 ns deadtime. 
We verify that the ratio of detected photon counts vs. total triggering count reaches 9.5%. 
At this ratio at which, based on coherent state detection probability calculation, the average 
photon number is \i — 1. The power meter and the id-201 were both employed for the initial 
calibration only, and these two independent approaches agree each other within the error 
bars in the DE measurement results. In subsequent measurements, we use the attenuator 
only to set the photon flux. 

The photon event, detected by the NFAD, is registered by a time tag unit whose time 
resolution is 156 ps. The time window width was set to be 1 ns for all of our measurements, 
unless explicitly mentioned. Each measurement was carried out for a second and repeated 
10 times. Fig. [5] shows the DE measurement results using Eq. [2} The experimental DE 
exhibits a saturation at darkcount rate above 10 kHz. The maximum detection efficiency 
for the NFADs is ~ 14%, and, with reasonable darkcount rates, it can reach <~ 10%. 

The saturation of detection efficiency (at about 14 %, illustrated in Fig. [5]) is actually due 
to the run-away increase in apparent DCR by afterpulsing. However, there is no way to sep- 
arate the random darkcount events from the afterpulses that they induce. The afterpulsing 
process is highly non-linear, since afterpulse events can induce additional afterpulsing events. 
Hence, there is always a point at which afterpulsing appears to increase quite dramatically 
with any additional bias. Nonetheless, it is expected that NFADs can provide significantly 
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FIG. 4. Measurement of the darkcount rates for the two NFAD devices. Both are measured at 193 
K with the same discriminator voltage settings. 
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FIG. 5. DE vs Dark count measurements results for two NFAD detectors. The discriminator 
threshold voltage level is set to 100 mV, and the testing temperature is 193 K. 

higher DE values (> 25%), similar to that of the id-201 detector, when operated in the 
gated mode with a sufficient hold-off time, although this imposes a reduction in detector 
"availability", that is undesirable for some experiments. 

Note that this detector efficiency is at the system level; no normalization was used to 
take into account the device active area, nor the insertion loss of the coupling path. In the 
DE test for both NFADs, we also vary the discriminator voltage levels. No significant DE 
improvements can be found by choosing different threshold values. 

We verify that the detectors act as bucket detectors. If this hypothesis is correct, then 
the detector click probability (the ratio of the detection event counts within the correlation 
time window to that of the total pulsed photon rate) is expected to behave as f{ji) = 



9 



T ! ! ' — II I 




0.01 0.1 1 10 100 0.01 0.1 1 10 100 

Average Photon Number \x Average Photon Number n 

FIG. 6. Verification of bucket detector model. The hollow squares are the experimental results. 
The solid lines are the simulation results based on modeling the detector as a bucket detector. 
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FIG. 7. Effect of temperature on detection efficiency and dark counts. The data points represent 
experimental measurements, and the lines are simple smooth fits. 

1 — exp(— r]v - n). The experimental results are compared with theoretical curve in Fig. [6j 
The curve is computed with the best fitting value of DE, ry D . The measurement results 
for two different NFADs are obtained at the darkcount rate approximately at 100 CPS. 
The other testing conditions are kept the same the DE measurement, except for the photon 
triggering rate, which is increased to 10 kHz. The good agreement for both detectors indicate 
that the bucket detector model is valid for our NFADs. This description is thus appropriate 
for modeling further quantum optics measurements. 

A key benefit of our unique readout circuit is that we operate the NFAD at lower temper- 
ature than most InGaAs SPADs are operated at (typically around 233 K). Fig. [7] shows the 
effect of temperature on the relation between DE and DCR. Once the DE and its associated 
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darkcount rate are measured, we can find the detector noise equivalent power (NEP) via 
NEP = ^V%D- The best noise equivalent power is estimated to be 8.1 x 10" 18 W Hz~ 1/2 
for Detector 1; and 1.4 x 10" 17 W Hz~ 1/2 for Detector 2 at 193 K. We will also study the 



figure of merit, which depends on the detector jitter time, in Section V 



B. DE measurements using a continuous wave pumped SPDC source 

Many quantum communication protocols are based on entangled photon pairs created 
from SPDC. It is of interest to characterize NFADs using photons from a SPDC source 
for several reasons: the non-Poissonian photon distribution of SPDC source behaves with a 
thermal-like property^, Most importantly, correlated photons allow to verify the detection 
efficiency with least assumption on photon flux calibratioii^, and avoid such calibration 
uncertainties in the WCP scheme; Moreover, the photon pairs are created at random times, 
allowing us to convincingly demonstrate the free-running mode of NFADs. 

The unique property of SPDC photon pairs has been exploited to characterize single 
photon detectors ever since the seminal work by Burnham et a/. 33 , other researchers have 
published their approache d 32 | 34 | 35 J to reduce the uncertainties in determining the absolute 
detector efficiencies. The experimental accuracy is limited by a wide varieties of nonideal 
factors^. The typical experimental configuration for characterizing a SPDC source operated 
with a continuous-wave (CW) laser, uses two free running single photon detectors, and time 
correlation photon counting logics. Note that the production times of photon pairs are 
random in this case. 

An alternative configuration!^, which is widely used particularly when the down converted 
photon pairs in the 1550 nm wavelength range, can be the use of of one free running detector, 
in combination with a gated SPAD, which typically has a higher efficiency. 

In our experiment, we tested two combinations of detectors, as displayed in Fig. [8| 
The SPDC source is an integrated waveguide based on a periodically poled lithium niobate 
(PPLN) substrates. Fig. [8] (a) illustrates the first setup, which is comprised of two free 
running NFADs, with their detection events recorded by a time tag unit. Both free-running 
NFADs have the figure of merit H > 1 x 10 7 at A = 1550 nm, as defined in Ref. [381 

Fig. [8] (b) shows the commonly used configuration: a free running NFAD is employed to 
trigger the second gated id-201 ULN SPD module. Its settings are: detection efficiency is 
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FIG. 8. Experimental setups used for two testing the NFAD 1 using SPDC. In dashed box (a), two 
free running NFAD detectors. In dashed box (b), one NFAD is placed in one of photon arms, a 
gated commercial detector (IDQ) is at the other arm. Photon detections from the NFAD are used 
to trigger an id-Quantique commercial single photon detector. 
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FIG. 9. (a) The inferred detection efficiency of the NFAD, from a direct comparison with an id-201 
commercial detector. The solid line (hollow circle) is the inferred DE using the SPDC scheme at 
1530 nm; the dash line (cross circle) is the inferred DE corrected for stray light detection; the solid 
line (square) is the DE characterized by the WCP scheme at 1550 nm. (b) The darkcount observed 
in the presence ( "apparent") and absence ("intrinsic") of stray light coupled into the optical fiber 
connecting to the detector when varying the bias voltage. 



10%, deadtime is 40 /xs, and gate width is 2.5 ns. The SPDC in PPLN is set to be in the 
non-degenerated mode, where the two down converted photons have different wavelength 
of (1570 ± 12.5) nm and (1530 ± 12.5) nm, respectively. The photons are separated by a 
dichroic mirror and fed into two detectors connected to channels 1 and 2. 

The overall procedure is based on two successive measurements depicted in Fig. [8] (a) 
and (b). For the two free running detector scheme, we keep NFAD 1 at 6.1 x 10 3 counts 
per second on average. To achieve this condition, we maintain a constant bias voltage to 
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NFAD 1 and a constant pumping power to the PPLN. We vary the bias voltage to the 
second detector NFAD 2. At each voltage step, we record the total single counts, as well as 
coincidental counts using the time tag unit. 

The net count (total single rate minus darkcount rate) of the down converted photons 
in channel one (CH-1) and channel two (CH-2) are given by N\ = rjiT]mN(l + N 2 = 
V2VD2N(l + k 2 )] their coincidental counts N c is given by N c = r]iT]mV < 2 r lD2N , where N is 
number of photon pairs, rji and rj 2 are the efficiencies of the two channels, and rjm and t]d 2 
are the detection efficiencies of the two detectors; Ki and k 2 are the afterpulsing coefficients 
of NFAD 1 and NFAD 2, which will be addressed in the following section. Therefore, from 

§ NFAD 
(a), we can find the ratio Xi = ^ = N c /N\ using two NFADs; 

Subsequently, if we switch to the scheme in Fig. M (b), we can find another ratio X2 = 



In practice, the afterpulsing coefficient is hard to maintain constant when the darkcount rate 



- N'JN[. Therefore, we can infer the efficiency of NFAD 2 via r^ AD = Vd^Xi/X2- 



is high. Additionally, the accuracy of this detector efficiency measurement is also limited by 
wide varieties of other nonideal factor^. 

Fig. [9] (a) displays the inferred detection efficiencies ^f AD of NFAD 2, and the results are 
compared with the aforementioned WCP method. The two plots are at the same bias voltage 
and discriminator threshold settings. The calculation of r]^ AD 18 based on the assumption 
that the channel 2 efficiency r] 2 , and NFAD 1 afterpulsing coefficient K4, are both kept 
constant in two consecutive schemes. The DE inferred from the SPDC method is plotted 
against apparent darkcounts when the stray light photons are included. Those stray light 
photons impose no effect on the inferred DE values, but reduce the SNR of the NFADs. We 
performed a separate darkcount measurement to estimate the number of darkcounts caused 
by stray light shown in Fig. [9] (b). This allows us to correct the inferred DE values based 
on the substraction of stray light detection. The corrected DE (the dash line in Fig. [9] (a)) 
exhibits improved SNR, and the overall shape of the WCP and SPDC methods agrees. 



IV. AFTERPULSING 

Afterpulsing is the most significant side effect for InGaAs SPADs; it introduces extra 
counts, in addition to the intrinsic photon and noise counts of the detector. It is mostly 
owing to defects in the semiconductor materials form carrier traps which hold the avalanche 
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charge carriers. The lifetime of the trapped charge carriers ranges from a few ns up to some 
tens of /is in conventional InGaAs/InP SPADs working in the gated Geiger mod^Sl. As the 
operating temperature is lowered, such as the case in this work, the afterpulsing problem 
becomes more prominent because the trapping is longer. 

The NFADs afterpulsing measurements are performed with the same setup as for the DE 
characterization (see Fig. 2). First, we measure the photon response using a fast digital 



oscilloscope instead of the time tagger. The measured screen shot is displayed in Fig. 10 (a). 
We also simulate the voltage across the intrinsic SPAD inside the NFAD with the SPICE 
circuit model (refer to Fig. [3] (b)). The transient simulation results is shown in Fig. [To] (b), 
where the time constant is determined to be 1.1 /is (R — 1.1 Mfi, (7 = 1 pF). 

We use the time tagger to measure the detector recovery. To identify which event is the 
photon detection, we measure a synchronization signal from the function generator on CH-1 
(see Fig. [2]) as the time reference. This enables us to discriminate between photon detections 
and events caused by dark counts or afterpulsing. All of the afterpulsing measurements are 
taken with the NFAD darkcount rate ^ set to about 100 CPS . Because we only record time 
tags for detector click event where the amplitude is higher than a preset threshold voltage, 
the weak avalanche events for the first few tens of ns after a detection are not recorded. 

Fig. [To] (c) illustrates this effect as an effective detector deadtime, following the rising 
voltage of NFAD recovery process right after a photon detection event. Depending on the 
operating temperature and the discriminator voltage settings, we observed that the NFAD 
effective deadtime is at least 40 to 50 ns long (dash line 1). After the deadtime period, we 
find a rapid recovery of afterpulsing events, rising within a time of about 25 ns (between 
dash line 1 and line 2). 



In Fig. 10 (d), we quantitatively measure the lifetime of the trapped avalanche carriers 
using a standard exponential decay function 41 , Coexp[— (t — Xd ea d)/7d]- r^ad, in our case, 
is the NFAD effective deadtime resulting from the finite threshold of the time tag unit. 
Theoretically, this value is around zero. With higher temperatures, the effective deadtime is 



slightly shorter, as shown in Fig. [TOJ is the characteristic detrapping time. Our best fit 
value to this curve is — 1.1 /is at both 193 K and 203 K, while the prefactor Co is smaller 
at 203 K. If we compare t& with voltage recovery time constant, it turns out that they are 
almost identical, which implies that detrapping of the residual avalanche charge is mainly 
induced by the detector recovery voltage. 
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FIG. 10. (a) Accumulated measurements of the transient detector response; this plot is the screen 
capture from a fast oscilloscope collected for 6 hours; the saturated amplitude is about 300 mV. 
(b) The simulation for the transient voltage drop across an intrinsic SPAD inside the NFAD. Both 
plots share the same time scale of 5 /xs. (c) Measurement of the detector deadtime and recovery 
time. The points are the measured photon detections in a 1 ns time resolution, and the solid line 
is the guide to eye. (d) Afterpulsing decay measurement and fitting. The random release of the 
trapped charges induced by photon detection events results in the afterpulsing probability decaying 
exponentially. The simulation curve is the solid line and the scattered dots are the accumulated 
measurement results for two temperatures. 



Fig. [IT] (a) shows the afterpulsing behavior of our NFADs for different discriminator volt- 
ages. We notice that increasing the threshold voltage can reduce the amount of afterpulsing 
for a short time after the detection, but the effect on the overall number of afterpulsing 



events is small. Fig. 11 (b) illustrates the effect of operating temperature on afterpuls- 



ing probability. As expected, afterpulsing becomes more prominent at lower temperatures. 
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FIG. 11. (a) Afterpulsing under two threshold levels of the time tag unit input. The temperature is 
at 193 K. (b) Afterpulsing probability density at three distinct temperatures, while the darkcount 
rate was kept constant at about 100 CPS. All data are collected for 1000 s. The threshold was set 
to be 0.1 V. The reference clock rate, used to synchronize the time frame, is set to 10 kHz. 
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FIG. 12. Afterpulsing parameter, ft, measurement at 193 K by WCP characterization setup for 
both two NFADs. (a) The afterpulsing coefficient as a function of the darkcount rate; (b) The 
afterpulsing coefficient as a function of the detection efficiency. The wavelength of probing WCP 
pulses is 1550 nm for all measurements. 



However, it is important to note that the measurements are performed at constant darkcount 
levels, and therefore the detector efficiency is lower for the higher temperature settings. 

Another interesting approach to model afterpulsing behavior of a NFAD is by introducing 
afterpulsing coefficient defined as k = Nt ~^ c ~ D , where N c is the correlated count rate, N t is 
the total raw count rate, and D is the darkcount rate at a given bias voltage. We vary the 
bias voltage to collect the raw counts and compute ft. The results are displayed in Fig. [12 
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(a) and (b). 

The comparison of the afterpulsing performance of two detectors indicates that Detector 
I has a better afterpulsing performance than that of Detector II. We attribute this to two 
factors: (1) Detector II has a larger active area than Detector I. Consequently, there are 
more defect trap sites in the active region that give rise to afterpulsing effects, by trapping 
carriers from one avalanche and releasing them at a later time when the detector becomes 
re-armed. (2) The larger area has a larger capacitance. This results in a proportionally 
larger number of charges generated when the capacitance is discharged with each avalanche. 
This increased charge flow per avalanche gives rise to more trapped carriers at defect sites, 
and thus the afterpulsing effect increases since more carriers are trapped and subsequently 
det rapped. The combination of these two effects creates a significant rise in afterpulsing 
events as the detector size is increased. The unwanted effects can be alleviated using active 
hold-off circuitry, which is widely used in gated mode detectors. However, we chose to keep 
the circuit as simple as possible, in order to avoid complexity of the electronics, its addition 
of stray capacitance, and undesirable elevate detector temperature. 

Note that the coefficient k is obtained via the WCP scheme. As we mentioned in the 
previous section, due to changes in experimental conditions including surround light sources, 
the observed darkcount rates in the SPDC experiments were different than for the WCP 
scheme. To estimate k using darkcount rate for the SPDC scheme, one must take into 
account the joint detection efficiency. 

V. TIME RESOLUTION 

For quantum optics applications, the timing jitter of the detector, r^ter, is another key 
parameter. The jitter time of the detector results from the fluctuations of its response 
time. Many mechanisms contribute to this variation. One source of jitter is the readout 
electronics. Another is intrinsic to the detector, primarily due to the stochastic build-up 
time of avalanche carriers within the multiplication region^ when the SPD detects a photon. 

Here, we measured the jitter time using a Becker & Hickl SPC-130 PCI card. We change 
the gain of the time to amplitude converter (TAC) to set the time resolution to 1.22 ps 
per division. This setting is used for all of the measurements. The other key technique for 
measuring high precision jitter time is to carefully adjust the zero crossing (ZC) level of both 
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inputs, usually set close the threshold levels. 



TABLE I. Timing jitter measurement!^ at 10 KHz repetition rate of pulse laser, 0.1 photon and 
1.0 photon on average per pulse. The measurements are accumulated for 1000 s. 



Low Limit 


ZC Level 


Total Count 


Peak Count 


FWHW 


Dark Count 


Average 


(mV) 


(mV) 






(ps) 


(CPS) 


Photon number 


-150 


-86 


100261 


799 


154 


~ 150 


0.1 


-150 


-96 


95620 


1186 


92 


~ 150 


0.1 


-140 


-96 


96172 


1163 


69 


~ 200 


0.1 


-110 


-96 


88066 


1361 


70 


~ 150 


0.1 


-100 


-96 


88118 


1405 


59 


~ 150 


0.1 


-90 


-96 










0.1 


-150 


-96 


706596 


15780 


31 


~ 150 


1.0 



The measurement setup is illustrated in Fig. [13] (a). The results are plotted in both 



logarithmic and linear scales by Fig. 13 (b) and (c), respectively. We took measurements 
for two different values of the average photon number per pulse. The low average photon 
number, 0.1 photon per pulse, effectively simulates the regime of single photon detection. At 
this level of /i, the jitter time is mostly below 150 ps, which is compatible to the resolution 
of our time tag unit. This photon flux level simulates the case present in many quantum 
optics experiments, such as detecting photon pairs generated by nonlinear optical processes. 
We also set \± = 1, not only a comparable photon flux to what is used for decoy state QKD 
protocol, but also to get an indication of the timing jitter intrinsic to the electronics. In 
all of the tests, we maintained the darkcount rate below 200 CPS. 

The measurement settings and jitter measurement are also summarized in Table [TJ For 
/i = 1, the measured jitter is 31 ps. For /i = 0.1, the timing jitter becomes wider also 
due to the width of laser puls^, which is specified to be <300 ps. Most measured jitters 
are around 100 ps, depending on the setting of PC-130 card ranging, ranging from 59 ps 
to 154 ps. Most importantly, all of the detection events fall within a 400 ps time length, 
regardless of measurement settings, which by itself is a very good timing response. The 
typical response of SPADs, including silicon ones, is that their time distribution has a long 
tail on the order of a few ns (see Ref. [38]) even if their FHWM-jitter is a few tens of ps. 
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FIG. 13. Measurement of the detector jitter by Becker & Hickl SPC-130 time-correlation card 
with a time resolution of 1.22 ps. (a) Experimental setup; (b) Measurements plots in logarithmic 
scale; (c) Normalized measurement plots in linear scale. The data are taken with 0.1 (sold lines) 
and 1.0 (dash line) photons per pulse for 1000 s, with different combinations of settings to the 
constant-fraction discriminator (CFD), the start signal for the time measurement. The function 
generator (AFG 3252), which produces the stop signal, is connected to the "Synch" input of the 
card. The detector output is connected to CFD input after an pulse inverter transformer. The 
average photon number per pulse, from 0.1 to 1, is controlled by changing the variable attenuator 
(VOA) settings. 
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Table [I] also suggests a trade off between the total counts and the combination of "Low 
limit" and "ZC level" settings, by which we can collect higher total counts, but at the cost of 
longer jitter time. However, no long recovery tail is found in all of our jitter measurements. 
Lastly, based on these jitter measurements, the detector figure of meritP^ can be calculated 
from the expression H = j^ 2 ^ — . Our observations show that the best figure of merit is 

■Lt'T jitter J 

6.3 x 10 7 for Detector 1. 



VI. SHOWING FEASIBILITY OF ENTANGLED PHOTON QKD OVER 
400 KM TELECOM FIBER 



Entanglement-based QKD can achieve the longest distance when the photon pair source 
is symmetrically put in the middle of two parties Alice and BofcPES Ma et a/. 46 first showed 
this configuration can tolerate up to 70 dB channel loss. Later Scheidl et alW^ performed the 
experiment and verified that the highest loss of their entanglement QKD achieved 71 dB. 
Here, we used a similar source, i.e. a CW pumped SPDC photon pairs. The quantum 
channel loss is simulated using a programable optical attenuator. The experimental setup is 



illustrated by Fig. [T4J To demonstrate the feasibility of a real QKD system, our setup has 
the same optical components as a standard QKD receiver, therefore we consider all potential 
losses from Alice, Bob and the entangled photon source. 



PPLN Waveguide 



OOP 



Toptica Laser 



PBS 



NFAD 1 Time Tag 




FIG. 14. Experimental setup for the investigation of High Loss QKD. Photon pairs are produced 
in a PPLN waveguide, and are sent through a variable attenuator. They are then split up by a 
dichroic mirror (DM) and detected by the NFADs. PBS: polarizing beam splitter; Attenuator: 
programmable optical attenuator; HWP: half wave plate; M: mirror; Col.: collimator. 

We performed 68 hours of measurements to test three attenuation values set by the 
variable optical attenuator. We use a time tag unit to collect the counts from two NFAD 
detectors. The time resolution can be chosen in the postprocessing stage. The internal 
optical losses are characterized for two wavelengths, averaging 8.22 dB per arm. The total 
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FIG. 15. The accumulated counting results collected by a 16-channel time tag unit with a time 
resolution of 156 ps, at three levels of attenuation settings, for (a) 62.4 dB and 24 hours collection 
time; (b) 66.4 dB and 20 hours collection time; (c) 70.4 dB and 24 hours collection time. Solid blue 
lines are Gaussian distribution fittings, with each of them having cr=0.25 ns , which is dominated by 
the timing resolution of the time tag unit, (d)-(f) are the experimental results (left) in comparison 
with the simulation results (right). In each of them, all above three levels of channel loss are taken 
into account, (d) accidentals-SNR: the ratio between the correlated and accidental counts; (e) the 
total correlated counts; (f) the total accidental background counts. 



loss is determined to be 62.4 dB, 66.4 dB and 70.4 dB by separate measurements. The 



darkcount rate for both detectors is set to be below 100 CPS. Fig. [15] displays the summary 
of our experimental results in (a)-(c), as well as measurements in comparison with to the 
simulated results in (d)-(f). In the simulation, the detector efficiencies are chosen to be 
Vdi = 8% at 1570 nm wavelength and rjm = 6% at 1530 nm wavelength. 
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The analysis of the experimental data and simulation of secret key rate is based on the 
framework of Ref. [471 However, we only focus on the accidentals visibility by looking at the 
ratio between the correlated and accidental counts. We also assume our source is an SPDC 
source with its characteristic emission statistics, operating at the average photon pair rate 
TV, is very low (4.5 MHz in our experimental measurements). t\a and t\b are the channel 
total efficiencies, each of them contains four contributions: (1) the finite detector efficiency 
(Vdi, Vd2)] (2) PPLN source coupling loss; (3) receiver optical loss (t]ca, Vcb) and (4) link 
attenuation loss (tsA, ^sb), which is set by the programable attenuator. So tja = VcAtsAVDi, 
and t]b = VcBtsBVD2- In general, tja <C 1 and t]b <C 1 due to photon transmission losses. 
Based on these conditions, the lower bound of quantum bit error rat€P^ Qber is- 

Qber = {2 + - [(i - — - - ^ - — VB + m } ■ (3) 

where the average photon pair number per time window \ip — NW. Here, W is the time 
window, and TV is the average photon pair rate. It is limited by the detector jitter time as 
well as the resolution of the time tag unit. From the previous experiments, we determined 
that jitter time of the NFADs is shorter than 150 ps, which is shorter than the time res- 
olution of the time tag unit 156 ps (from DotFast / UQDevices^) . Note that the Becker 
& Hickl card, which has much better time resolution, could not be used in this applica- 
tion because it is not suitable for the continuous time stamping of photon pair detections 
over the long accumulation times we require (> 20 hours). Consequently, in all of our 
measurements and simulation, the smallest time window is chosen to be 156 ps. Based on 
this timing resolution, TV could be as high as 1 GHz and still provide a 9:1 ratio for the 
1-photon pair probability versus the 2-photon pair probability. We assume ka = k>b = 2 



as the afterpulsing coefficients to fit both detectors (measurements shown in Fig. 12), given 
the detector efficiencies in conjunction with their darkcount settings. Da and Db are the 
darkcount rates, and are both set to 100 CPS. 

Using the QBER obtained from Eq. [3j we can estimate the secret key rate, i?, by: 

R > \ {N VA rj B [l-f {Q BER ) H 2 (Qber) - H 2 (Qber)}} , (4) 

where H 2 is the binary entropy function, H 2 (x) = — xlog 2 (x) — (1 — x) log 2 (1 — x). The 
correction factor is / (Qber) = 1-22. 
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FIG. 16. The estimated secret key rate for entanglement QKD versus the total channel loss. 
In the numerical simulation: for the cases of k = 2, we use detector efficiency of 8% and 6%, 
respectively; and for the cases of n — 1, two detector efficiencies are chosen to be 4% and 3%, 
respectively. The solid squares are based on actual measurements of the signal to noise for a given 
loss. The lines show secret key rates calculated for different source rates and afterpulsing coefficient. 
Two additional x-axes show corresponding QKD distances according to the total channel loss, for 
standard SMF-28 single mode fiber and ultra low loss fiber, SMF-28 ULL. These results indicate 
that entanglement-based QKD over 400 km is feasible. 



We experimentally simulated secret key rates for two types of telecom optical fibers: 
a standard telecom fiber (SMF-28) whose loss, a^, is 0.2 dB per km, and an ultra low 
loss (ULL) telecom fiber (SMF-28-ULL) 50 , whose loss is 0.164 dB per km. Note that in 
both cases, we have introduced 6 dB extra loss, to account for optical losses including 
interconnections at Alice, Bob and central source sites. The results are showed in Fig. [THJ 
The discrete dots in the figure are based on the experimental data for three levels of optical 
attenuations. The results suggest that the distance of entangled pair based QKD is feasible 
over 300 km via the standard telecom fiber, and can reach over 400 km in telecom ULL 
fiber. One caveat is that the finite key size effect has not been taken into account for both 
simulations. Nonetheless, given the time resolution of our detector and the possibility to 
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achieve an entanglement photon source whose rate can reach up to a few hundreds of MHz, 
entanglement-based QKD system over 400 km is possible. 

Another remarkable observation is that normally a detector free of afterpulsing (k = 0) 
will have its cut-off QKD distance limit proportional to the photon pair rate. This is not 
the case for the current NFAD design, because our detector is free running and no extra 
afterpulse suppressing circuit has been introduced, leading to its afterpulsing coefficient of 
k — 2. With this nonzero value, the QKD cut-off distance is no longer proportional to N. 
Instead the secret key rate vanishes when > 70 MHz approximately. However, one remedy 
to this situation is to decrease the current detector efficiency to 4% and 3%, respectively. 
As a result, the afterpulsing parameter drops to less than one. Consequently, in the case of 
k — 1, the secret key rate does exist even when N = 100 MHz. However, as indicated in 



Fig. 16, its secret key rate is lower than the case when N = 20 MHz, and the efficiency of 
both detectors are twice as high with k = 2. Nevertheless, the regime where k = 1 will be 
beneficial for studying novel high rate SPDC entanglement photon pair source. 



VII. CONCLUSION 

We have reported a novel design to read out single photon responses for NFAD oper- 
ating at low temperatures. The system circuit diagram has been presented in detail. We 
carried out a SPICE model simulation to compare the conventional readout method with 
our approach. The simulation result agrees with the experimental measurement. Detection 
efficiencies of both NFADs are characterized using TCSPC measurements under different 
temperatures and threshold settings via WCP source. Subsequently, we have used two free 
running NFADs and compared them with an id-201 detector. The inferred DE of NFAD II 
is derived from coincidence counts of the id-201 matches the expected NFAD efficiencies 
obtained in WCP characterization. By scanning the average photon number in WCP, we 
confirmed that the NFAD acts as a bucket detector^ i n the quantum optics model. We also 
found that the detection repetition rate can be close to 100 kHz without DE degradation. 
The afterpulsing was tested by varying the threshold settings and ambient temperatures. 
The temporal distribution of afterpulsing events follows an exponentially decaying trend. 
Our measurement of the NFAD jitter time can achieve a FWHM of 31 ps. The dark noise 
below the rate of lO^ns" 1 for the best NEP 8.1 x 10~ 18 W • Hz" 1 / 2 . The fi eure of merit for 
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SPDs was determined to reach 6.3 x 10 7 . We also demonstrated that an entanglement-based 
QKD system can tolerate higher losses up to 70 dB, which infers successful key distribution 
distance over 400 km of ULL telecom optical fiber. 

We like to stress that the operating temperature for NFAD is maintained in continuously 
cryogenic free mode with a simple electronic circuitry. No cryogenic facility is involved; 
no vacuum is needed. The cooling power of our system (10 W) allows it to be an readily 
expandable platform to include higher numbers of NFADs, where the heat emitted from each 
detector is less than 1 mW, even for count rates of about 1 MHz detection pulse rate per 
detector. These advantages enable this NFAD system to be a good candidate in quantum 
optics, quantum cryptography, and SPDC entanglement source characterization applications 
for the telecom band. 
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